Quantum Speed Limit for Perfect State Transfer in One Dimension 
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The basic idea of spin chain engineering for perfect quantum state transfer (QST) is to find a 
set of couphng constants in the Hamiltonian, such that a particular state initially encoded on one 
site will evolve freely to the opposite site without any dynamical controls. The minimal possible 
evolution time represents a speed limit for QST. We prove that the optimal solution is the one 
simulating the precession of a spin in a static magnetic field. We also argue that, at least for solid- 
state systems where interactions are local, it is more realistic to characterize the computation power 
by the couplings than the initial energy. 
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The transfer of a quantum state from one part of a 
physical unit e.g. qubit, to another part is a crucial in- 
gredient for many quantum information processing pro- 
tocols Q. In fact, the ability to transmit "flying qubits" 
is known as one of the DiVincenzo's desiderata ^] for 
quantum communication. For long distance communica- 
tion, it is natural to rely on optical means However, 
for systems involved in a sufficiently small scale where 
qubits can interact either directly or indirectly through 
other qubits, it seems to be more natural to exploit the 
interactions between them directly. In these scenarios, a 
typical assumption usually made is that the interactions 
between pairs of qubits can be arbitrarily switched on 
and off [J, . In this way, a quantum state can be trans- 
mitted by a series of swap (SOS) operations. However, 
even in the absence of decoherence, the intrinsic problem 
of this protocol is that each swap would necessarily intro- 
duce some errors due to dynamical controls, the fidelity 
of the transferred states would therefore decay exponen- 
tially as the number of swaps increases. It is therefore of 
great practical interest to study schemes of perfect state 
transfer by free Hamiltonian evolution which is primar- 
ily designed to minimize dynamical controls as much as 
possible. 

This work is also motivated by the question of the ulti- 
mate speed limit of solid-state quantum computers. The 
Margolus-Levitin theorem suggests that the shortest 
time for a quantum state to evolve into an orthogonal 
state is limited by the initial energy. With this result, 
one can quantify the speed limit for elementary opera- 
tions involving one or two qubits However, this the- 
orem does not give an appropriate bound on the speed 
limit for the initial state to evolve into a particular state, 
for example the one dimensional quantum state transfer 
(QST) problem [cf. Eq. 0]. The aim of this paper is to 
establish the speed limit for the QST problem. Following 
the scheme described in Ref . , one can generalize the 
study to more general quantum operations. In contrast 
with the Margolus-Levitin theorem, the quantum speed 
limit for QST will be characterized directly by the cou- 
plings in the Hamiltonian, instead of the initial energy. 



As we shall see, this makes it not only experimentally 
more accessible, not also physically more reasonable. 

Imperfect QST |^, 0| over a uniform spin chain was 
studied numerically, showing that perfect state transfer 
is impossible for long chains. However, it has been sug- 
gested that perfect QST, without dynamical con- 
trols, is possible if we allow the couplings between the 
spins to be non-uniform. Moreover, it has also been 
pointed out |lOl | that there are many possible solutions 
for constructing spin chains which allow perfect state 
transfer, and the searching for the solutions is in gen- 
eral an inverse eigenvalue problem. Recently, there have 
been many quantum information processing proposals 
[Il[IlIIlIlE[IllIl[iS.I9,2a,i23 related to the con- 
cept of spin chain engineering. To boost the performance 
of the spin chains, or to evaluate the ultimate speed limit 
of these proposals, it is highly desirable to know how to 
engineer spin chains with optimal performance. 

In the following, we shall show that the optimal solu- 
tion for perfect state transfer is exactly the one proposed 
in Ref. . Amusingly, this optimal set of couplings 

is the same as that describing spin precession in a static 
magnetic field. In discussing the one dimensional QST 
problem, one usually considers the Hamiltonian for the 
eng ineered spin-1/2 chains to be of the following form 



N-l 

H=y^ 



where crj's are the standard Pauli matrices for the site j, 
and N is the total number of sites in the spin chain. 
Both and Aj are real constants to be determined. 
We shall adopt the convention that |0j) (|lj)) represents 
spin-down (spin-up) state. 

In the quantum state transfer protocol, one site is en- 
coded with a quantum state a |1) |0) representing the 
information of one qubit, while the rest are initialized to 
be spin-down. Define U (t) = e^*^"^, where r is the evo- 
lution time and Since U (t) |000...00) = |000...00), 
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the state is said to be perfectly transferred whenever 

C/(t) I100...00) = e'"^ I000...01) , (2) 

where is some known phase. Although the trans- 
ferred state in general differs from the input state by 
a relative phase factor, i.e. ae*"^ |1) + /3|0), the phase 
factor could be corrected by some local operations and 
is thus usually ignored. To investigate further, it is 
clear that we can simply focus on the single-particle 
subspace {|1) = |100...0) , |2) = |010...0) , ...}, where |a;), 
X = 1,2, 3,.., TV refers to a state with a single excitation 
at site X. In this subspace, the Hamiltonian in Eq. |^ is 
real, symmetrical and tridiagonal: 



Hs = 



/ Ai wi 

OJl A2 0J2 

U)2 A3 

V 







ljn~i Aat 



(3) 



The eigenvalues, which must be non-degenerate [25j, and 
the corresponding eigenvectors are denoted by Ek and 
jcfc) respectively. Let us consider not only transferring 
states from one end to the other but more generally the 
transition amplitude between the state \ x) and its mirror- 
inverted state \x), with x=N—x+l, 



N 



{x\UiT)\x) =J2{S^\ek){ek\x)e-'''''^ 



(4) 



fc=i 



Provided that the Hamiltonian in Eq. or Eq. (jSJ is 
mirror symmetrical, i.e. A^ = Aj and ujj = luj^^j, one 

can show 0| that {x\ek) = (— 1)'^ (ccjefc), and hence 
|(S|;7(r)|x)| = lif 



= (-1) e 



fc 10 



(5) 



The phase factor e*"^ here is exactly the same as that 
appeared in Eq. ||2Jl. It has been shown [l3| that the 
condition in Eq. ||SJ) is not only sufficient but also neces- 
sary for perfect state transfer in mirror symmetric chains. 
Given an eigenvalue spectrum {Ek\ satisfying Eq. 
the corresponding coupling constants {ujj , Aj } can be 
uniquely determined. The task of solving for the solu- 
tions of \ujj , Xj } is therefore an inverse eigenvalue prob- 
lem |2J,|25j , which generally has to be solved numerically. 
A spin chain is said to be engineered if its couplings are 
found in this way. 

Perhaps the simplest solution to the Eq. ^ is the 
linear spectrum |l2l |. By mapping the states |a;) to 
be the eigenstates of the Jz angular momentum oper- 
ator, this solution also describes the precession of a spin 
J—{N—l)/2 under a constant magnetic field pointing 
along the x-direction. Explicitly, the matrix elements 
of Hs can be chosen as 



A, = and 



1 



VjW- 



(6) 



TABLE I: Comparing the efficiency 77 (defined in Eq. ||7|l) 
for various possible solutions to the problem of perfect state 
transfer through a spin chain of spins. For convenience, 
we have scaled the maximum tunneling matrix element tJmax 
such that the minimal evolution time r = tt is fixed. 



Type of Spectrum 



Efficiency rj 



Linear" 
Quadratic'' 
Gapped linear L 
Gapped linear IL* 
Cosine*^ 



^max 

q {2a+N) Cj^ 
(l+4g/Ar)J;n 

g 

108.5 



1 

q (2a+N) 
(l+4g/iV) 

q 

14.0 



1. (imax is defined in Eq. |2J and = 1 by definition. 
[. For iV^>l, q (o+l) >l/2 as required by _E;i-_Eo>1, 
positive integer and o is a rational number, 
j. For N'^ > 1. 
I . For N = i only. 

|. For A'' = 31 (no analytic solution given). 



"See Ref.n: 
''See Ref.E 
where g is a 



<^See Ref. 
''See Ref. 
■^See Ref. 



The eigenvalue spectrum is linear in the sense that the 
energy spacing is uniform, Ek= — {N—\)/2-\-k—\, k — 
1, 2, 3, TV. From Eq. ^ and (O, perfect state transfer 
can be achieved for a time period of r = tt. 

Apparently, there are unlimited number of eigenvalue 
spectra satisfying the condition in Eq. ((SJ). A number of 
particular solutions have also been found 0, 0, 0, 
recently. To compare different sets of inter-qubit coupling 
constants {ujj} generated from different energy eigen- 
value spectra {Ek}, the efficiency rj can be quantified by 
the evolution time r for completing the task of perfect 
state transfer, subject to the constraint that the maxi- 
mum value of inter-spin coupling Wmax = maxjwj} be- 
ing normalized. Alternatively, one may fix the evolution 
time and compare the maximum coupling strength for 
different chains. Both cases can be properly captured by 
defining the efficiency as 



UJn 



(7) 



where Wmax and f are respectively the maximum coupling 
and the evolution time of a reference spin chain, which 
we shall choose to be the one described in Eq. 10, i.e. 
f = TT and 



UJ„ 



1 for odd iV, 
for even iV. 



(8) 



With these definitions, our goal is to show that 77 > 1 for 
all engineered spin chains. This implies that the set of 
couplings in Eq. ij^l is optimal and the quantum speed 
limit for one dimensional QST is established (i.e. given 
the same maximum couplings Wmax = ^-max, t = rfT > 
f). A comparison for different schemes of perfect state 
transfer is given in Table 

We start with the observation pointed out in Ref. . 
Suppose the eigenvalues are ordered &s Ei > E2 > ■ ■ ■ > 
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i^TV- Define Ak = Ek ~ Ek+i and the range of the spec- 
trum as Ae = El- En ^ 'I2k=i ■ From Eq. ©, the 
evolution time is limited by the minimum energy inter- 
val Amin = min{Afe}, r > 7r/Ai„i„, while f = 7r/Amin- 
From Eq. Q, we have 



> 



An 



-^max^min 



(9) 



^max^min 



Suppose we scale the eigenvalue spectra such that A,„i„ — 
Amin, then ?7 > 1 if Wmax > <^inax- Wc shall show that it 
is indeed true. 

In the single-spin subspace, applying the Hellmann- 
Feynman theorem to the Hamiltonian in Eq. we 
have dEk/dXj — ^ (ek\ [a^ + l) le^) and dEk/dujj ~ 
2(efe|i)(i+l|efe). It is clear that dEk/dXj > 0. On the 
other hand, it is known ^] that the eigenvector |efe) has 
exactly fc — 1 sign changes in the position basis In 
other words, dEi/dujj > and dEN/duj < 0. Suppose 



we replace all ujj 



and A, 



Amax, An 



respec- 



tively, and invoke the solution of the eigenvalues for a 
uniformly coupled chain: A + 2ti> cos(-j^^), we have 



Ae < Aa +4tj 



(10) 



where A^ = Amax — Amin is the range of the spatial varia- 
tions of the local potentials. Suppose, in the absence of, 
or with uniform, local potentials, we set A^ = 0. Under 
the constraint that the evolution time to be the same, i.e. 
Amin > Ainin = 1, the range of any nonlinear spectrum 
must be at least greater then that of the linear spectrum 
by one, Ae > Ae + I ~ N , a.s the ratio between any two 
Afc's must be a rational number One can then show 
from Eq. ((HJ and (fllHl that ujmax > Wmax- Of course, this 
argument is based on the assumption that the terms Xj 
play no role in the eigenvalue spectrum except a constant 
shift. The question is: could the speed hmit be boosted, 
if spatially varying local potentials are allowed? The aim 
of the following paragraphs is to exclude such possibility, 
and at the same time provide a more rigorous proof. 

Exploiting the symmetry of our problem, we can di- 
vide the Hs in Eq. © into two different subspaces 
For even N, they are simply + and 

75 - = l,2,..,iV/2. The X iV Hamiltonian 

reduces into two 



N N 

X ones: 



( Xi uoi 

UJi A2 

L02 A3 



\ ^N/2^^N/2 J 



(11) 



These two matrices are almost identical except one 
diagonal matrix element Xfq/2^^N/2- Recall that 
{x\e].) = (—1)''' {x\ek), the eigenvectors of Hs are 



automatically grouped into these two subspaces. In 
other words, the eigenvalue spectra of these two ma- 
trices are respectively = {Ei,E^,Ez, ...,E]si-i\ and 
/ifc — {E2, E4, Eq, E]\j}. Consequently, the trace 
difference between the two matrices gives C1J7V/2 = 
(l/2)I]fe (^fe ^ Aife) = (l/2)I]fc ^2fe-i, which is mini- 



mized when all A2fc-i 
the case ujf^/2 



which corresponds to 



'^max- Hence we can conclude that 

^max > > a),nax for CVCU N. 

For odd iV>5 (the case N—3 can be analytically 
solved separately), it becomes slightly more compli- 
cated. The subspaces are spanned by — 

and + together with \m), where j — 

1, 2, 3, m— 1 and m = (A^-fl) /2. The matrix formed 
by the latter is 






V 











\ 



Am-2 
Wm-l A„j_i \/2ll!„ 
V2uj,n Xm 



(12) 



and the counterpart by the former can be obtained 
from the above matrix by removing the last row and 
the last column. Similar to the case of even N, the 
eigenvalue spectra are i>k — {Ei, E^, Ej^} and /ife = 
{E2, i?4, Ei^^i} respectively. Our goal is still to show 
'-^m ^ '^max whcn Amin ^ Amin- To procccd, ouc Can 
diagonalize the upper part of the matrix in Eq. I|12|l . i.e. 
the part from Ai to A„i_i , and the resulting matrix is of 
the bordered diagonal form (arrowhead matrix) 



( Xi ••■ 
A2 



b2 



A,„_i bm-i 
\bi 62 ■ ■ • byn-i a J 



(13) 



where a — Xm — J2T=i — X^fcLi A^fc the trace differ- 
ence, and bk are the off-diagonal matrix elements satis- 
fying the condition 2w^ = Y^^=i b^, as a result of of the 
diagonalization. 

The characteristic polynomial P (E) — JlfcLi " ^fe) 
of the matrix above can be written as 



P [E) = {E-a)Q (E) - 



in—l 

E 



bl 



-{ {E-t^k) 



Q{E) , (14) 



where Q {E) 
have 



]Xk=i {E - Hk)- By setting E = fik, "^e 



U"lik if^k - Mj) 



> 



(15) 
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where the negative sign here is because of the interlac- 
ing property Vk > tJ-k- One can solve for Um recursively 
using Eq. H14|l and (|15|l . We found a recursive relation 

is the counterpart of LOm when we solve the same in- 
verse eigenvalue problem without the last two eigenvalues 
i?7V-i and En. Since uj'^ is a positive sum of the prod- 
ucts AjAk, it is minimized when all = Amin, and 
hence cj,nax > '-^m > Wniax- This completes our proof for 
rj > 1 for all possible engineered spin chains. 

In the following, we shall analyze the implications of 
the speed limit. It is natural to ask, apart from the 
advantage of not requiring dynamical controls, if the 
schemes of engineered chains can be fundamentally faster 
than the series of swaps (SOS) protocol mentioned at 
the beginning. The answer is positive. Suppose we now 
consider the ideal case where each swap operation can 
be achieved with the same maximum possible coupling 
f^max in Eq. ©, the evolution time for each swap is 
''"0 ~ 7''/2(Di„ax- Since it will totally take N — 1 steps 
to transfer a state through a chain of N spins, we can 
generalize our definition of the efficiency for this case as: 
Vsos = — 1) /2ujnia-^, which is an increasing function of 
N, and one can easily show that 1 < ?7sos < 2 for N > 2. 
Therefore, the engineered spin chains can have roughly a 
factor of two gain in speed for sufficiently long chains. 

According to the Margolus-Levitin theorem 0|, the 
minimum time for a quantum state ji/'i) to evolve into an 
orthogonal state \ipf), where {1pf\^pi)=0, is limited by the 

initial energy: max ( 2(5-^0) ' 2Ke) ' where E={'>p,\H\Tp,) 
is the initial energy, Eq is the ground state energy and 
AE~y/ {tpi\{H — is the energy uncertainty of the 

initial state. The computational power of physical sys- 
tems as quantum computers was analyzed |17| along this 
line of thought. Here we shall argue that, at least for 
solid state systems where the interactions are local, it is 
more realistic to characterize the computational power 
by the couplings in the Hamiltonian than the initial en- 
ergy. We shall again focus on the ID QST problem [cf. 
Eq-@] as an example, as the generalization to more com- 
plicated quantum operations is straightforward through 
the scheme proposed in Ref. For the linear spec- 

trum, we have E—Eq=N/2 for large N . Given the same 
amount of initial energy N/2, the minimal evolution time 
for the SOS protocol is t = tt/N (as E-Eq = AE). For 
a linear chain of N spins, the SOS protocol seems to 
be as efficient as the engineered chains. However, as we 
have just analyzed, for each swap the tunneling strength 
would be N/2, which is twice of the maximum possible 
coupling a)niax assumed before! 

In conclusion, we have generalized the study of quan- 
tum speed limit to the problem of quantum state transfer 
in one dimension. We have argued that for solid state sys- 
tems it is more realistic to characterize the computation 
power explicitly by the coupling terms in the Hamiltonian 



than the initial energy. This work also implies that quan- 
tum algorithms, when implemented in solid state devices, 
can be optimized, if we employ the concept of spin chain 
engineering More practically, systematic errors are 
expected to be significantly minimized with the reduc- 
tion of dynamical controls. We lastly comment that such 
speed limit, if not confined with nearest-neighbor inter- 
actions in one dimension, could be overcome by increas- 
ing either the network complexity or the Hamiltonian 
complexity Details to be discussed elsewhere. 
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